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Abstract. We have calculated the rate of the decay 6 — > dVi to second order in in the 
limit that the h and c quarks have equal masses. The results here confirm recent calculations 
done in the opposite limit where the c-quark is much lighter than the 6-quark. 



1. Introduction 

In the effort to find physics beyond the Standard Model, one of the methods being implemented 
is over-constraining the unitarity triangles that arise from the Cabbibo-Kobayashi-Maskawa 
(CKM) quark mixing matrix. One important quantity that must be measured for this is the 
value of the CKM parameter \Vcb\- This is typically done in one of two ways: inclusive or exclusive 
measurements. In exclusive measurements specific decay rates, for example B D*~^i~17i, are 
measured and compared with theoretical predictions to extract a value for \Vcb\- The inclusive 
measurements, however, look at all decays from a meson with a 6-quark to a daughter meson 
containing a c-quark. Currently, the most accurate measurement we have for \Vcb\ is from 
the inclusive decays. The agreement between the two measurements, however, is only about 
2.5%. As such, one would like to improve on these measurements to aid in over-constraining the 
unitarity triangles. 

Looking specifically at the inclusive decays, the current value of \Vcb\ is. 



with a large portion of the error coming from higher orders in the perturbative calculations 
[1]. This determination of \Vcb\ included the leading and next-to- leading order perturbative 
calculations along with an order of magnitude estimate of the second order corrections called 
the Brodsky-Lepage-Mackenzie (BLM) correction. 

Naturally, the next step is to calculate the full second order corrections. This was done in two 
separate calculations: one done numerically [2] and one analytically [3]. The analytic calculation 
involved an expansion in the 6 to c quark-mass ratio, p = rric/mb, with a light c-quark. This 
expansion gave results that were consistent with the numerical results of when the ratio of 
masses was in the region of the physical value. If one wanted to go above this value though, 
the analytic result started to diverge. This is solely due to the truncation of the expansion in 
p that was done. In order to extend our knowledge of the second order corrections to /3 = 1 we 
have carried out the same calculation but in the opposite limit of the mass ratio i.e. p ~ 1 [4j. 
This expansion provides an independent check of the results in [3] as well as providing accurate 
values for the decay rate through the full range of the quark-mass ratio. 
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Figure 1. A sample of diagrams that needed to be calculated for the second order corrections 
to the decay width. Dashed lines denote leptons and curly lines denote gluons, while thick and 
thin lines denote b and c quarks, respectively. 

2. Methods 

To calculate the second order corrections to the decay width of 5 — > ciUi, we treat the masses 
of the c and b quarks to be almost equal and expand in small parameter 5 = 1 — p. The second 
order corrections involve calculating decays with up to five daughter particles. We found it 
convenient to use the optical theorem and instead calculate the imaginary parts of four loop 
diagrams, some of which are shown in Figured! One of the benefits of calculating loop diagrams 
is that we can then use all of the machinery that has been developed. In particular, for this 
problem, the method of asymptotic expansion simplifies the calculation immensely [5jj [6] . 

Asymptotic expansion is a method of separating the large and small momentum contributions, 
or regions, of the loops in a diagram so that the resulting integrals are more tractable. In 
this calculation, we happen to be in a kinematic configuration that was ideal for asymptotic 
expansion. When applied, the region with all loop momenta large ends up not contributing to 
the decay rate. This region typically gives the most difficult integrals and so our calculation 
was greatly simplified. Another benefit of asymptotic expansion in this calculation was that all 
diagrams could be reduced to five different integrals with at most two loops. 

The other important property of loop integrals used, is that of recurrence relations. This is 
a method based on integration-by-parts identities [71 18] that gives relations between Feynman 
integrals and it is often the case that the resulting integrals are easier to calculate than the 
original one being considered. Specifically, we used the Laporta algorithm [9l [TO] to treat the 
diagrams that have a three gluon vertex: e.g. Figure [Ija. We used two programs, ROWS [11] 
and FIRE [12] in order to do this reduction. 

With these simplifications in place, we were able to obtain an expansion that is highly 
convergent and surprisingly easy to compute, especially compared to the expansion from the 
opposite limit that was done in [3]: see also [TO] . 

3. Results 

The total width is customarily written in the form 



with Gp the Fermi constant, and m;, the pole mass. We have Cp = 4/3, Ca = 3, Tp = 1/2 in 
QCD and n/ = 3 the number of light flavours in our theory. Full analytic results for Xq and Xi 
can be found in [14J, while analytic expansions for X2 can be found in [SjH]- 

The final result of our calculation is an expansion in 6 with Xc, Xb and Xi terms calculated 
up to 0{5^^), Xp to 0(5^'^), and Xa to 0{5^^). As an example of the results, we provide the 





X2 = CpXp + CaXa + TpimXi + Xc + Xb) 



(3) 



Xfe term to order 5"^ here. This term receives contributions from diagrams hke that in Figure [2]3. 
The full results can be found in the source files for the preprint of 



/184 32 2\ ,5 / 8 2\ .6 /107444 3848 .7 

Note that the decay is highly suppressed as it starts at 5^ . As well, the first two terms include 
fractions that are simpler than the 6"^ term. This occurs in all of the second order terms and has 
an interesting explanation as to why. These two terms come from the zero recoil limit where 
rric = mi) and the mass of the leptons is zero. This calculation has been done in [15] and we 
have verified the first two terms in our expansion using these results. 

If we look at the plot of as compared to the expression in [3] , Figure [21 we can see the 5^ 
suppression. Figure [2] also shows how well our expansion agrees with the previous one and we 




Figure 2. Plot of the contribution as a function of p = 1 — 5. The solid line is our expansion, 
while the dashed line is the expansion from [3]. 



were able to reproduce the minimum at about p = 0.2. 

A plot of the full result is shown in Figure [3j In this plot, we can see that the two expansions 
agree well in the region p k, Q.l — 0.4. As well, our expansion is accurate to about 2.5% all the 
way to /9 = 0, (5 = 1, where our expansion isn't expected to work well. Given the simplicity of 
this calculation and this fact, it is conceivable that a calculation of the 0{a^g) corrections could 
be computed with results that extend across the full range of mass ratio rric/mf,. 



4. Conclusions 

We have calculated the second order corrections to the decay b ciui in the limit of a heavy 
daughter quark. This expansion agrees with a previous expansion in the opposite limit [3j and 
verifies some interesting properties of the previous result. As well, our result is highly convergent 
and thus gives accurate values for the decay rate at any value of the expansion parameter 

5. Along with the simplification of the calculation through use of procedures like asymptotic 
expansion and the Laporta algorithm, this gives promise for higher order calculations in the 
future. 




Figure 3. Plot of the full second order correction as a function of p = 1 — (5. Again, the solid 
line is our expansion, while the dashed line is the expansion from [3j. 
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